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Abstract 



In this paper we investigate the Wightman function, the renormalized vacuum expecta- 
^v^j \ tion values of the field square, and the energy-momentum tensor for a massive scalar field 

. with general curvature coupling inside and outside of a cylindrical shell in the generalized 

spacetime of straight cosmic string. For the general case of Robin boundary condition, by 
using the generalized Abel-Plana formula, the vacuum expectation values are presented in 
the form of the sum of boundary-free and boundary-induced parts. The asymptotic behavior 
of the vacuum expectation values of the field square, energy density and stresses are inves- 
tigated in various limiting cases. The generalization of the results to the exterior region is 
Oh. given for a general cylindrically symmetric static model of the string core with finite support. 



PACS numbers: 03.70.+k, 98.80.Cq, 11.27.+d 



1 Introduction 

Cosmic strings generically arise within the framework of grand unified theories and could be 
produced in the early universe as a result of symmetry braking phase transitions ^ Al- 
though the recent observational data on the cosmic microwave background radiation have ruled 
out cosmic strings as the primary source for primordial density perturbations, they are still 
candidates for the generation of a number interesting physical effects such as the generation of 
gravitational waves, high energy cosmic strings, and gamma ray bursts. Recently the cosmic 
strings attract a renewed interest partly because a variant of their formation mechanism is pro- 
posed in the framework of brane inflation [HllUOi]- In the simplest theoretical model describing 
the infinite straight cosmic string the spacetime is locally flat except on the string where it has a 
delta shaped curvature tensor. In quantum field theory the corresponding non-trivial topology 
leads to non-zero vacuum expectation values for physical observables. Explicit calculations have 
been done for various fields [H]-j2Hl- The case of quantum fields at non-zero temperature was 
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also investigated [2U1 [23 EM EU • Vacuum polarization effects by the cosmic string carying a 
magnetic flux are considered in Refs. |321 1331 Another type of vacuum polarization arises 
when boundaries are present. The imposed boundary conditions on quantum fields alter the 
zero-point fluctuations spectrum and result in additional shifts in the vacuum expectation val- 
ues of physical quantities, such as the energy density and stresses. In particular, vacuum forces 
arise acting on constraining boundaries. This is the well-known Casimir effect (for a review see, 
|35j ) . In this paper we will study the configuration with both types of sources for the vacuum 
polarization, namely, a cylindrical boundary coaxial with a cosmic string assuming that on the 
bounding surface the field obeys Robin boundary condition. For a massive scalar field with an 
arbitrary curvature coupling parameter we evaluate the Wightman function and the vacuum 
expectation values of the field square and the energy-momentum tensor in both interior and 
exterior regions. In addition to describing the physical structure of a quantum field at a given 
point, the energy-momentum tensor acts as the source in the Einstein equations and therefore 
plays an important role in modelling a self-consistent dynamics involving the gravitational field. 
The vacuum densities for a Robin cylindrical boundary in the Minkowski background are in- 
vestigated in [HE] (see also [SZI for the case of two coaxial cylindrical surfaces). In Ref. [2^] a 
cylindrical boundary with Dirichlet boundary condition is introduced in the bulk of the cosmic 
string as an intermediate stage for the calculation of the ground state energy of a massive scalar 
field in (2+l)-dimensions. 

We have organized the paper as follows. The next section is devoted to the evaluation 
of the Wightman function for a massive scalar field in a generalized cosmic string spacetime 
in both interior and exterior regions of a coaxial cylindrical boundary. In section [21 we present 
simple formulae for the vacuum expectation values of the field square and the energy-momentum 
tensor in the corresponding boundary-free geometry. By using the formula for the Wightman 
function, in section [I] we evaluate the vacuum expectation values of the field square and the 
energy-momentum tensor inside a cylindrical boundary. Various limiting cases are investigated. 
In section El we consider the corresponding quantities for the outside region. The generalization 
of the results for the exterior region in the case of a general cylindrically symmetric static model 
of the string core with finite support is given in sectional Finally, the results are summarized 
and discussed in section [7| 

2 Wightman function 

2.1 Bulk and boundary geometries 

In this paper we consider a scalar field ip propagating on the {D + l)-dimensional background 
spacetime with a conical-type singularity described by the line-element 

N 



with the cylindrical coordinates (x 1 , x 2 , . . . , x D ) = (r, (f>, Z\ , . . . , zn), where iV = D — 2, r ^ 0, 
^ 4> ^ —°o < z i < +°° and the spatial points (r, c/>, zi, . . . , zn) and (r, <p + 4>o, Z\, ■ ■ ■ , z^) 
are to be identified. In the standard D = 3 cosmic string case the planar angle deficit is related to 
the mass per unit length of the string /i by 2ir — cfto = 87rG/i, where G is the Newton gravitational 
constant. It is interesting to note that the effective metric produced in superfluid 3 He — A by a 
radial disgyration is described by the D = 3 line element (|2,1|) with the negative angle deficit, 
that is (fio > 2tt, which corresponds to the negative mass of the topological object [3*5] . 

For a free massive field with curvature coupling parameter £ the field equation has the form 






(2.2) 
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where Vj is the covariant derivative operator and R is the scalar curvature for the background 
spacetime. The values of the curvature coupling parameter £ = and £ = £d = (D — 1)/4D 
correspond to the most important special cases of minimally and conformally coupled scalars, 
respectively. We assume that the field obeys Robin boundary condition on the cylindrical surface 
with radius a, coaxial with the string: 

A + B—)<p = 0, r = a. (2.3) 



dr 

Of course, all results in what follows will depend on the ratio of the coefficients in this boundary 
condition. However, to keep the transition to Dirichlet and Neumann cases transparent, we 
use the form 1)2. 3 j) . In this section we are interested in the corresponding positive frequency 
Wightman function in the regions inside and outside of the cylindrical surface due to the fact 
that the vacuum expectation values (VEVs) of the field square and the energy-momentum tensor 
are expressed in terms of this function. In addition, the response of a particle detector in an 
arbitrary state of motion is determined by this function (see, for instance, [SH])- By expanding 
the field operator in terms of a complete set of eigenfunctions {ip Q (x), <p* a (x)} satisfying the 
boundary condition and using the standard commutation relations, the Wightman function is 
presented as the mode-sum 

{0\p(xMx')\0) = M*)<P*a&), (2-4) 

a 

where a is a collective notation for the quantum numbers specifying the solution and |0) is the 
amplitude for the corresponding vacuum state. The form of the eigenfunctions is different for 
the regions inside and outside the cylindrical shell and we will consider these regions separately. 

2.2 Wightman function in the region inside the shell 

In the region inside the cylindrical surface the eigenfunctions satisfying the periodicity condition 
on 4> = 4>Q are specified by the set of quantum numbers a = (n, 7, k), n = 0,±1,±2, 
k = (/bi,..., kjy), —00 < kj < 00, and have the form 

(p a (x) = f3 a J q \ n \{^r) exp [iqncj) + ikry - iut) , (2.5) 

u = vV + k 2 + m 2 , g = 2vr/</)o, (2.6) 

where rn = (zi, . . . ,2jv) and Ji(z) is the Bessel function. The eigenvalues for the quantum 
number 7 are quantized by the boundary condition 1)2.3)1 on the cylindrical surface r = a. From 
this condition it follows that for a given n the possible values of 7 are determined by the relation 

7 = W». j = l,2,---, (2.7) 

where X n j are the positive zeros of the function «7 g [ n i(z), J q \ n \{\ n ,j) = 0, arranged in ascending 

order, X n j < A n j-fi, n = 0, 1,2, Here and in what follows, for a given function f(z), we use 

the notation 

f(z) = Af(z) + (B/a)zf'(z). (2.8) 

It is well-known that for real A and B all zeros of the function J q \ n \{z) are simple and real, 
except the case Aa/B < —q\n\ when there are two purely imaginary zeros. In the following we 
will assume the values of Aa/B for which all zeros are real. 

The coefficient (5 a in 1)2.5)1 is determined from the normalization condition based on the 
standard Klein-Gordon scalar product with the integration over the region inside the cylindrical 
surface and is equal to 

ol _ ^n,jT qn (X n j) 

Pa (2vr)^ a 2 ' { ] 
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where we have introduced the notation 



T qn (z) 



1 2 n 2 ) J 2 qH {z) + z*J>^{z) 



(2.10) 



Substituting the eigenfunctions ()2.5|) into the mode-sum formula ()2.4|) with a set of quantum 
numbers a = (rejk), for the positive frequency Wightman function one finds 



(oi^M^lo) 



^ke ik(r " 



II y, cos[qn{(j) — 4>)] 

n=0 



oo 



x > \PlJ q n{ir)JUir')e-^ t ~ t \=x n , 3 /a, 



(2.11) 



where the prime means that the summand with n = should be taken with the weight 1/2. As 
we do not know the explicit expressions for the eigenvalues X n j as functions on n and j, and 
the summands in the series over j are strongly oscillating functions for large values of j, this 
formula is not convenient for the further evaluation of the VEVs of the field square and the 
energy-momentum tensor. In addition, the expression on the right of (|2.11j) is divergent in the 
coincidence limit and some renormalization procedure is needed to extract finite result for the 
VEVs of the field square and the energy-momentum tensor. To obtain an alternative form for 
the Wightman function we will apply to the sum over j a variant of the generalized Abel-Plana 
summation formula (301 



if 00 IF 



dz ■ 



Iqn(z) 

e- qn7Ti f(zef)+e gnwi f(ze-f) 



(2.12) 



where h{z) and Ki{z) are the modified Bessel functions. This formula is valid for functions f(z) 
analytic in the right half-plane of the complex variable z = x + iy and satisfying the condition 
|/(z)| < e(x)e Cl \ y \, c\ < 2, for \z\ — > oo and the condition f(z) = o(z 2q ^ n ^ 1 ) for z — > 0, where 
e(x) — > for x — > oo. By taking in formula (|2.12() qn = 1/2 and B = 0, as a special case, we 
obtain the Abel-Plana formula. 

To evaluate the sum over j in (|2.11|) as a function f(z) we choose 



zJgn(zr/a)J qn (zr'/a) 
a/ A; 2 + m 2 + z 2 /a? 



cxp 



-iy/k 2 + m 2 + z 2 /a 2 (t - t') 



(2.13) 



Using the asymptotic formulae of the Bessel functions for large values of the argument when 
n is fixed (see, e.g., |3B)> we can see that for the function f(z) given in ()2.13j) . the condition 
to formula ()2.12j) to be satisfied is r + r' + \t — t'\ < 2a. In particular, this is the case in the 
coincidence limit t = t' for the region under consideration, r, r' < a. Formula H2.12|) allows to 
present the Wightman function in the form 



(0\<p{x)<p(x')\0) = (0 s \<p{x)<p(x')\0 s ) + (<p(x)<f(x')) a , 



(2.14) 



where 



(0 s \<p(x)cp{x')\0 s 



I f d N k Mru-r',,) 



<Ao7 (2tt) 



N 



dz 



ze 



-i(t-t')Vz 2 +k 2 +m 2 



V z 2 + k 2 + m 2 



y~] cos[qn((j) -4>)\ J qn (zr)J qn (zr'), 



(2.15) 



n=0 
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and 



d N k ik(r|| _ r;) "*> ^ cosh 



(t - t')V^ 2 - k 2 



m 



ncpo J (2tt) n A/F+w 2 ' V-z 2 - fc 2 



2 



CO 



x V'cosfen^ - ^)]J.p.W^(^) ^ft (2-16) 



n=0 



(za) 



In the limit a — > oo for fixed r,r', the term ((p(x)ip(x')) a vanishes and, hence, the term 
{Q s \ip(x)(p(x')\Q s ) is the Wightman function for the geometry of a cosmic string without a cylin- 
drical boundary with the corresponding vacuum state |0 S ). Consequently, the term (ip(x)(p(x')) a 
is induced by the presence of the cylindrical boundary. Hence, the application of the generalized 
Abel-Plana formula allows us to extract from the Wightman function the part due to the string 
without a cylindrical boundary. For points away from the cylindrical surface the additional part 
induced by this surface, formula (|2.16|) . is finite in the coincidence limit and the renormalization 
is needed only for the part coming from ()2.15|) . 

2.3 Wightman function in the exterior region 

Now we turn to the region outside the cylindrical shell, r > a. In general, the coefficients in 
Robin boundary condition (|2.3f) for this region can be different from those for the interior region. 
However, in order to not complicate the formulae we use the same notations. The corresponding 
eigenfunctions satisfying boundary conditions (|2.3j) are obtained from Q2.5|) by the replacement 

Jqn{ir) -> g qn (>yr,ja) = J qn {ir)Y qn (-fa) - J qn (^a)Y qn (^r), (2.17) 

where Y qn (z) is the Neumann function. Now the spectrum for the quantum number 7 is contin- 
uous. To determine the corresponding normalization coefficient (3 a , we note that as the normal- 
ization integral diverges in the limit 7 = 7', the main contribution into the integral over radial 
coordinate comes from the large values of r when the Bessel functions can be replaced by their 
asymptotics for large arguments. The resulting integral is elementary and for the normalization 
coefficient in the region r > a one finds 

Pi = S 2i:)2 ~ D \ (2.18) 

20o u [ J\\ (7a) + Y q \ (7a)] 

Substituting the corresponding eigenfunctions into the mode-sum formula ()2.4|) . the positive 
frequency Whightman function can be written in the form 

i(? - *V7 2 + k 2 



„ 1°° d ^9qn {ir, ja)g qn ( 7 r', 7a) ex P 



J 2 n (~fa)+Y 2 n ( ia ) V7 2 + k 2 



(2.19) 



To find the part in the Wightman function induced by the presence of the cylindrical shell we 
subtract from (|2.19|) the corresponding function for the geometry of a cosmic string without a 
cylindrical shell, given by Q2.15[) . In order to evaluate the corresponding difference we use the 
relation 

^T^TT' - *M W) -~t, V-^MH^), (2.20) 
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where Hqn(z), s = 1, 2 are the Hankel functions. This allows to present the Wightman function 
in the form (|2.14|) with the cylindrical shell induced part given by 



(ip(x)(p(x')) a 



1 

Wo 



(2VT) 



oo 



AT 



cos[qn((f> - (j)')] ^ 



n=0 



s=l 



<i 77 



x ffg( 7 r)^( 7 r / ) J gw ( 7 a) 
77^ + ^ ffW( 7 a) 



i{t' -t)^^ 2 + k 2 



(2.21) 



On the complex plane 7 , we rotate the integration contour on the right of this formula by the 
angle tt/2 for s = 1 and by the angle —n/2 for s = 2. The integrals over the segments (Q,ik) 
and (0, — ik) cancel out and after introducing the modified Bessel functions we obtain 



(ip(x)(f{x')) a 



d N k 



7T0O J {27r) N 



ik(ri| — r 



Jk 



oo z cosh 
dz 



(t - t')Vz 2 - k 2 



y~] cos[qn((j) - (j)')]K qn (zr)K qn {zr 



n=0 



V 7 ? 

K qn (za) 



k 2 



(2.22) 



As we see the boundary induced part of the Wightman function for the exterior region is obtained 
from the corresponding part in the interior region by the replacements I ^ K. Note that in 
formulae ()2.15|) , ()2.16|) and ()2.22|) the integration over the directions of the vector k can be done 
using the formula 



d N k 
(2tt) n ' 

for a given function F(k). 



I 



ikr 



F(k) 



dk k N - 1 F(k) JN/2 ~ 1 ^ 



(2tt) n / 2 



(ifclxl)^/ 2 - 



(2.23) 



3 Vacuum expectation values for a string without a cylindrical 
boundary 

In this section we consider the geometry of a string without a cylindrical boundary. The one- 
loop quantum effects of the scalar field in this geometry have been considered in a large number 
of papers. The VEV of the energy-momentum tensor for a conformally coupled D = 3 massless 
scalar field was evaluated in Ref. [Sj- The case of an arbitrary curvature coupling is considered 
in Refs. [51 1101 IT2*] . The integral representation for the Green's function for a massive scalar field 
is considered in Refs. [51 1171 l2"Tl I22j . The corresponding local zeta functions are discussed in 
Refs. jSniEEl- However, to our knowledge, no closed formulae have been given in the literature 
for the VEVs of the field square and the energy-momentum tensor in the case of a massive field 
and in an arbitrary number of dimensions. The mass corrections in the limit mr <C 1 have 
been considered in Refs. |221 125j . Below we will derive simple exact formulae for the VEVs of 
both field square and the energy-momentum tensor assuming that the parameter q is an integer. 
These formulae will also shed light on the qualitative behavior of the VEVs in the case of non- 
integer q. The Wightman function for the geometry of a string without cylindrical boundary is 
given by formula (|2.15|) . For a field with the mass m and in the case of N parallel dimensions 
along the string we denote this function as G^ N ^(x, x' , m). The following recurrence relation is 
a simple consequence of formula (|2.15j) : 

G^(x,x',m) = - dk N cos [k N (z N - z' N )} G^ N _^(x, x' , \J k 2 N + m 2 ). (3.1) 
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Of course, this formula takes place for any other two-point function (see Ref. |21| for the case of 
Green function). After the integration over the angular part of the wave vector k using formula 
(|2.23|) . the integral over k is evaluated on the base of formulae given in Ref. [32] ■ Assuming 
that |rii —rn| > \t — t'\, one obtains the following expression 



v l-v 



{O s \(p(x)<p(x')\Q s ) = V ' — / dz z (z 2 + m 2 ) Kv-^vy/? + m 2 ) 

2 V v -k v $o Jo 



oo 

X 

n=0 



cos[qn((f> - <j)')]J qn {zr)J qn (zr'), (3.2) 



|(r H -rj,) 2 - (t-0 2 | 1/2 - (3-3) 



where we have introduced the notations 

D - 1 

v = ~2~' v ~ ~ L \\ 
The corresponding formula in the case |ry — rji | < \t—t'\ is obtained from Q3.2|) by the replacement 

AVi -> (vr/2)e^^- 1 / 2 )^ 1) 1 . In special cases, formula Q3.2JI agrees with the results form Refs. 
|121 117j . The analog formula for the Green's function in the case of more general background 
of the spinning cone and for a twisted scalar field is given in |22| . For integer values of the 
parameter q the expression for the Wightman function can be further simplified by using the 
formula |42| (see also Ref. jllj 1 ) 

oo ^ q-l 

y~] cos[qn((j) - (p')]J qn (zr)J qn (zr') = — Jo (zui) , (3.4) 

n=0 ^ 1=0 

where 

= r 2 + r /2 — 2rr' cos(/0o + <j) — (/)'). (3-5) 

By using ()3.4|) in Eq. ()3.2|) and evaluating the integral over z with the help of a formula from 
|42j . we arrive at the following expression 



D-l 

m 2 



(0.1^)^)10.) = — . l V 2 ,,/ 2 - (3-6) 
(2tt) 2 i=Q (-u 2 + w 2 ) 

where v and v are defined by formulae ()3.3|) . The / = term in formula (|3.6|) coincides with the 
Wightman function of the Minkowski space without boundaries, and the Wightman function 
for the geometry of the string is a sum of q images of the Minkowski spacetime functions. We 
could obtain this result directly by using the method of images (see also |12l for D = 3 and 
D = 2 massless cases). 

The VEV of the field square is formally given by taking the coincidence limit x' — > x of the 
Wightman function. However, this procedure leads to a divergent result. In order to obtain a 
finite and well defined result, we apply the renormalization procedure which corresponds to the 
subtraction from the Wightman function in the geometry of the cosmic string the corresponding 
function in the background Minkowski spacetime. So the renormalized VEV of the field square 
reads 

^ 2 )[± = (0s|^|0 s )-(0 m |^ 2 |0m) 
m ' 1 K v (2mryi) 

(2^ hi < 2mr ^ ' 

yi = sin(7rZ/g). (3.8) 



where we used the notation 
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From (|3.7I) it follows that the renormalized VEV of the field square is positive everywhere. In 
the limit mr ^> 1/ sin(7r/g), the main contribution to the VEV (|3.6j) comes from the terms with 
I = 1 and I = q — 1 and to the leading order one finds 

„ , x m / 2 - 1 e -2mrsin(n/q) 

(to 2 \^ « — - (3 q^i 

\r /ren r A \n/9 ■ n /9 1 i \ ' V° - "^ 

with exponentially suppressed VEV. 

In the case of a massless scalar field we obtain from (|3.7[) the following result 

2 (s) r((£>-l)/2) ^ r-p 

)U = 77—5+1-—— 2^ y i • ( 3 - 10 ) 

(4tt)— r^" 1 i=1 

For the case D = 2 this agree with the result in Ref. |16j . From (|3 . 7|) it follows that in the case 
of a massive field the expression on the right of this formula determines the leading term in the 
asymptotic expansion of the VEV of the field square near the string. For odd values of D the 
summation on the right of formula (|3.1U[) can be done by making use of the formulae 

-r / \ ^n{x) + n 2 T n (x) , 2 ■ -2/ \ ■ -2 /n i i \ 

T n+2 [x) = , 12(1) = ? sin (gx)-sm x, (3.11) 

for the sum 

T n (x) =^sin- n (x + /vr/g). (3.12) 
i=i 

In particular, one has J 2 (0) = (<? 2 - l)/3 and in the case D = 3 we obtain the well-known result 
for the VEV of the field square in the massless case [HI For D = 5, by using recurrence 
relation 1)3.11(1 for the evaluation of Ti(x), we find li(0) = (q 2 — l)(q 2 + ll)/45 and 

y>W = (qa - 1) ^ + 11) ,g = 5. (3-13) 

\V /ren 2880ttV 4 ' V ' 

It is worth call attention to the fact that the resulting expressions are analytic functions of q 
and by analytic continuation they are valid for all values of q. 

Now we turn to the VEVs of the energy-momentum tensor in the cosmic string spacetime 
without boundaries. These VEVs can be evaluated by making use of the formula 



(0\T ik (x)\0) = lim ViV k (0\<p(x)<p(a/)\Q) + 



<(V(x)|0>, (3.14) 



and the expressions for the Wightman function (j3.fi j) and for VEV of the field square. For the 
non-zero components one obtains 



0\(«) 



m D+i «4 ( 2 K v+2 (2mryi) 



ren D+ 



tf>S = ^E(^?-l)%^, (3-16) 



(2t) » S" Pmn/l)' 



(2vr 



2 z=i 



K u+1 (2mryi) K u+2 (2mryi) 



(2mryi) v+l (2mryi) 1 



(3.17) 



S 



and for the components in the directions parallel to the string we have (no summation over i) 

(T;)£ n = (T°)£ n ,i = 3,...,D. (3.18) 

It can be explicitly checked that for a conformally coupled massless scalar field this tensor is 
traceless. Of course, as the background spacetime is locally flat, we could expect that the trace 
anomaly is absent. From the continuity equation V^T^ = for the energy-momentum tensor 
one has the following relation for the radial and azimuthal components: 

d r (rTl) = Tl (3.19) 

It can be easily checked that the VEVs given by (|3.1fi|) . (|3.17|) satisfy this equation. For ^ ^ 1/4 
the radial stress is negative and one has the relation (T| )rcn > — ^(Ti^ren- In particular, this is 
the case for minimally and conformally coupled scalar fields. In the limit mr 3> l/sin(7r/g) to 
the leading order we have the following relations 

(Zff)W „ (l-^m^in 2 ^)^ 2 )^), (3.20) 
(T|)W « m 2 [l-4£sin 2 (7r/g)]^ 2 >«, (3.21) 

where the asymptotic expression for the VEV of the field square is given by formula (|3.9ft . 
For the radial stress one has (T^ren ~ — (2~2)ren/ [2mr sin(7r/g)]. In the same limit the energy 
density is positive for both minimally and conformally coupled fields. In figure^ the vacuum 
energy density is plotted for massive D = 3 minimally and conformally coupled scalar fields, as 
a function of mr for various values of the parameter q. As it has been mentioned above, for 
large values mr the energy density tends to zero, being positive for both minimal and conformal 
couplings. For the latter case the energy density is negative near the string for q > 1 and, hence, 
it has a maximum for some intermediate value of mr. For a minimally coupled scalar the energy 
density near the string is negative for q 2 > 19 and positive for 1 < q 2 < 19. 




1 1.5 2 2.5 3 3.5 1 2 3 4 5 



Figure 1: Vacuum energy density divided by m D+1 , (Tq )ren/w' D+1 , for minimally (left panel) 
and conformally (right panel) coupled D = 3 scalar fields as a function on mr for various values 
of the parameter q (the numbers near the curves). 

The formulae for the components of the vacuum energy-momentum tensor in the case of a 
massless scalar field are obtained in the limit m — > (no summation over i): 

1{S) _ r(QD + l)/2) «=i /(Oft,,) 
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where 

fW(Vi) = (1 - 4£)(D - - 1, = 4^ 2 - 1, (3.23) 

and the vacuum stresses obey the relation (T 2 )^ = —D(Tl)i&. In particular, for the case 
D = 3, by using the expressions for 2^(0) and 2^(0) given before, we obtain the result derived 
in I12| . In the case of a massive scalar field, the expression on the right of Eq. (|3.22j) 
determines the leading term in the asymptotic expansion of the vacuum energy-momentum 
tensor near the string. 



4 Field square and the energy-momentum tensor inside a cylin- 
drical shell 

4.1 VEV of the field square 

We now turn to the geometry of a string with additional cylindrical boundary of radius a. Taking 
the coincidence limit x' — > x in formula (|2.14|) for the Wightman function and integrating over 
k with the help of the formula 

,w, f°° k s g(z)dz ir N / 2 „ ( N + s 1\ f°° , , 9 ?x , s , , 

d k/ yy ' — = — — —B[ — — ,- / dz(z 2 -m 2 ) 2 g(z), (4.1 

JVk^ v^ 2 -k*-m> T(iV/2) \ 2 '2)J m K ' 

where B(x,y) is the Euler beta function, the VEV of the field square is presented as a sum of 
two terms 

((V|0> = <0 s |p 2 |0 fl > + (p 2 > a . (4.2) 

The second term on the right of this formula is induced by the cylindrical boundary and is given 
by the formula 



^ = -¥ eT dzz ^ - ™ 2 ) ^ ^rrW. ( 4 - 3 ) 



n=0 

where we have introduced the notation 



7T 2 



For the points away from the cylindrical surface, r < a, the integral in l|4.3|) is exponentially 
convergent in the upper limit and the boundary-induced part in the VEV of the field square is 
finite. In particular this part is negative for Dirichlet scalar and is positive for Neumann scalar. 
Near the string, r <C a, the main contribution to (f 2 ) a comes from the summand with n = 
and one has 

H 

As the boundary-free renormalized VEV diverges on the string, we conclude from here that 
near the string the main contribution to the VEV of the field square comes from this part. In 
particular, on the base of the results from the previous section we see that for integer values q 
the VEV of the field square is positive near the string. 

The part (f 2 ) a diverges on the cylindrical surface r = a. Near this surface the main contri- 
bution into (|4.3|) comes from large values of n. Introducing a new integration variable z — ► nqz, 



(^) a «-_£i_/ dzz(z 2 -m 2 a 2 )^^. (4.5) 
2a u i <6 n ./„,„ In (z) 
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replacing the modified Bessel functions by their uniform asymptotic expansions for large values 
of the order (see, for instance, and expanding over a — r, up to the leading order, one finds 



( i P 2 )a ~ — — — -rrn — / dz — = } n u z e 



(q/2a) D -H25 B0 -l) f°° ^ z D ^ ^ ^-^-^(l-r/qjyT+F 

(25 B0 - i)r (gg) 
(47t) 2 (a — r) 



(4.6) 



This leading behavior is the same as that for a cylindrical surface of radius a in the Minkowski 
spacetime. As the boundary-free part is finite at r = a, near the boundary the total renormalized 
VEV of the field square is dominated by the boundary-induced part and is negative for Dirichlet 
scalar. Combining this with the estimation for the region near the string, we come to the 
conclusion that in this case the VEV of the field square vanishes for some intermediate value of 
r. 

Now we turn to the investigation of the boundary-induced VEV given by ()4.3j) . in the lim- 
iting cases of the parameter q. Firstly consider the limit when the parameter q is large which 
corresponds to small values of 4>o and, hence, to a large planar angle deficit. In this limit the 
order of the modified Bessel functions for the terms with n ^ in ()4.3[) is large and we can 
replace these functions by their uniform asymptotic expansions. On the base of these expansions 
it can be seen that to the leading order the contribution of the terms with n ^ is suppressed 
by the factor q^ D ~ l ^ 2 {r / a) 2q and the main contribution to the VEV of the field square comes 
from the n = term: 

^"-^ I™ dzz^-rn^tS^-llizr), q»l, (4.7) 
20o Jm Io(za) 

with the linear dependence on q. In the same limit the boundary-free part in the VEV of the field 
square behaves as q D ~ l and, hence, its contribution dominates in comparison with the boundary- 
induced part. In the opposite limit when q — > 0, the series over n in Eq. (|4.3j) diverges and, hence, 
for small values of q the main contribution comes from large values n. In this case, to the leading 

_/oo 

order, we can replace the summation by the integration: Q f(Q n ) ~^ (V?) Jo°° dxf(x). As 

a consequence, we obtain that in the limit q — > the boundary-induced VEV in the field square 
tends to a finite limiting value: 

( ( P)a = --^- dx dzz(z -m) - I x (zr). (4.8) 
2vr Jo J m lx(za) 

Now we consider the limiting case obtained when 4>q — > 0, r, a — > oo, assuming that a — r 
and acpQ = ao are fixed. This corresponds to the geometry of a single boundary in the spacetime 
with topology i?^ -1 ' 1 ) x S 1 and with ao being the length for the compactified dimension. We 
introduce rectangular coordinates (x n , x' 2 , . . . , x' D ) = (x, y, z\, . . . , zn) with the relations x = 
a — r, y = a<j) in the limit under consideration. In this limit, from the quantities corresponding 
to the geometry of a string without a cylindrical surface we obtain the vacuum densities in 
the spacetime i?^ -1,1 ) x S 1 . These quantities are well-investigated in literature and in what 
follows we will consider the additional part induced by the presence of the boundary at x = 
0. The corresponding vacuum expectation values are obtained from the expectation values 
{■■■)a- For this we note that in the limit under consideration one has q = 2ir/(j)Q — > oo, 
and the order of the modified Bessel functions in formula (|4.3|) for n ^ tends to infinity. 
Introducing a new integration variable z — > qnz, we can replace these functions by their uniform 
asymptotic expansions for large values of the order. In the term with n = the arguments 
of the modified Bessel functions are large and we replace these functions by the corresponding 
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asymptotic expressions. After straightforward calculations, the vacuum expectation value of the 
field square is presented in the form 

(0|^|0) = (^)(°) + (^)«, (4.9) 

where (v? 2 )^ is the vacuum expectation value for the topology R( 1,D ~ 1 ^ x S 1 without boundaries, 
and the term 



OC 



u n =0 m " 

with m n = \J 'm 2 + (2irn/ao) 2 , is induced by the presence of the plate at x = 0. In (|4.10j) 
the terms with n ^ correspond to the contribution of Kaluza-Klein modes related to the 
compactification of the y direction. The boundary induced VEV for the field square in the 
spaces with topology i?^" 1,1 ) x E with an arbitrary internal space S is obtained in [IS] as a 



limiting case of the corresponding braneworld geometry. It can be checked that formula (|4.1 
is a special case of this formula for E = S . 

4.2 VEV of the energy-momentum tensor 

The VEV of the energy-momentum tensor for the situation when the cylindrical boundary is 
present is written in the form 

(0|T tt |0) = (0 s |T ifc |0 s ) + (T ik ) a , (4.11) 

where {Tik) a is induced by the cylindrical boundary. This term is obtained from the corre- 
sponding part in the Wightman function, (ip(x)<f(x')) a , acting by the appropriate differential 
operator and taking the coincidence limit [see formula (|3.14|) ] . For the points away from the 
cylindrical surface this limit gives a finite result. For the corresponding components of the 
energy-momentum tensor one obtains (no summation over i) 

Y0 n _g Jm -lqn\Z(l) 



with the notations 



f rz (y) + [i + q ^f)f z (y) 



2 2 2 

y z - m z r z 2 



F^[f(y)] = 

F$[f(y)] = \f\y) + jf(y)f\y)-l(i + ^f)f 2 (y), (4.u) 
F$lf(y)] = ( 2 ^)[^+( 1 + ^) 



,2^2 

f(y) ' 



+ ^f\y) - -/(y)/ / (y)j(4.i5) 
y z y 



and 



F$[f(y)]=F$[f( y )},i = 3,...,D. (4.16) 

It can be checked that the expectation values Q4.12JI satisfy equation (|3. 19|> and, hence, the 
continuity equation for the energy- momentum tensor. The boundary- induced part in the VEV 
of the energy-momentum tensor given by Eq. (|4.12|) is finite everywhere except at the points on 
the boundary and at the points on the string in the case q < 1. As we will see below, unlike to 
the surface divergences, the divergences on the string are integrable. 

In the case q > 1, near the string, r — > 0, the main contribution to the boundary part 1)4. 12j) 
comes from the summand with n = and one has 

A D 



f°° dzz 3 (z 2 - m 2 a 2 ) ^ (*) (4.17) 

J ma Io{ z ) 
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with the notations 



= 2i - \ + Iz^/f!, F (0 (z) = e - I, i = 1, 2. (4.18) 

For q < 1 the main contribution to the boundary-induced part for the points near the string 
comes from n = 1 term and in the leading order one has 

297rr J (g + l) J m I ff («o) 



where 



F! (0) = g(2£ - 1/2), F[ 1] = if } = (2g - (4.20) 



As we see, in this case the VEVs for the energy-momentum tensor diverge on the string. This 
divergence is integrable. In particular, the corresponding contribution to the energy in the region 
near the string is finite. 

As in the case of the field square, in the limit g> 1 the contribution of the terms with 
to the VEV of the energy-momentum tensor is suppressed by the factor q( D ~ 1 ^ 2 (r/a) 2q and the 
main contribution comes from the n = term with the linear dependence on g. In the same 
limit, the boundary-free part in the VEV of the energy-momentum tensor behaves as q D+1 and, 
hence, the total energy-momentum tensor is dominated by this part. In the opposite limit when 
q <C 1, by the way similar to that used before for the VEV of the field square, it can be seen that 
the boundary-induced part in the vacuum energy-momentum tensor tends to a finite limiting 
value which is obtained from Q4.12JI replacing the summation over n by the integration. The 
described behavior of the VEVs as functions of q is clearly seen in figure|2]where the dependence 
of the boundary-induced vacuum energy density and stresses on the parameter q is plotted in the 
case of D = 3 minimally and conformally coupled massless scalar fields with Dirichlet boundary 
condition and for r /a = 0.5. 



o.i 




Figure 2: Boundary-induced parts in the components of the energy-momentum tensor, 
a D+1 (T-) a , evaluated at r/a = 0.5 as functions of the parameter q for minimally (full curves) 
and conformally (dashed curves) coupled D = 3 massless scalar fields with Dirichlet boundary 
condition. The numbers near the curves correspond to the value of the index i. The azimuthal 
stress for the case of conformally coupled field is obtained from the zero-trace condition. 

The boundary part (Zf) diverges on the cylindrical surface r = a. Introducing a new 
integration variable z — > nqz and taking into account that near the surface r = a the main 
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contribution comes from large values of n, we can replace the modified Bessel functions by their 
uniform asymptotic expansions for large values of the order. To the leading order this gives 

D(t - Zo)(25 B0 - 1) T (D + 1\ Z _ Q2 D (421) 

This leading divergence does not depend on the parameter q and coincides with the corresponding 
one for a cylindrical surface of radius a in the Minkowski bulk. For the radial component to 
the leading order one has (T±) a ~ (a — r)~ D . In particular, for a minimally coupled scalar field 
the corresponding energy density is negative for Dirichlet boundary condition and is positive for 
non-Dirichlet boundary conditions. For a conformally coupled scalar the leading term vanishes 
and it is necessary to keep the next term in the corresponding asymptotic expansion. As the 
boundary-free part in the VEV of the energy-momentum tensor is finite on the cylindrical 
surface, for the points near the boundary the vacuum energy-momentum tensor is dominated 
by the boundary- induced part. 

In the limit 4>q — > 0, r, a —* oo, with fixed values for a — r and cufio = do, proceeding in a 
similar way to that used for the field square, the VEV of the energy-momentum tensor can be 
written as 

(0|^|0) = (^)(°) + (^)«, (4.22) 

where (Tf)^ is the corresponding quantity for the topology w^' D ~^ x S 1 without boundaries, 
and the term (no summation over {) 

OO 

W>™ = £> - ^% )e ^, ,4.3) 

is induced by the boundary located at x = 0. In (|4.23j) . m n is defined by the relation given in 
the paragraph after formula (|4.1U|) and we have introduced the notations 

F?\z) = (4£-l), 2 + ^i (4.24) 

27m x 2 



F^(z) = (4e-l)^— j , (4.25) 

with F^\z) = 0. In particular, in this limit the boundary-induced vacuum stress in the direction 
perpendicular to the plate vanishes. Formula (|4.23|) is a special case of the more general result 
for the spaces with topology i?^ -1 ' 1 ) x E with an arbitrary internal space S obtained in 01] as 
a limiting case of the corresponding braneworld geometry. 

In figure |31 we present the graphs for the boundary induced parts of the components of the 
energy- momentum tensor as functions of r/a. The graphs are plotted for D = 3 massless scalar 
field with Dirichlet boundary condition on the cylindrical surface and for the cosmic string bulk 
with q = 4. In this special case for a minimally coupled scalar field, the boundary-induced 
vacuum energy density vanishes on the cosmic string. 

5 VEVs in the region outside a cylindrical shell 

In this section we consider the VEVs induced by the cylindrical boundary in the exterior region 
r > a. Taking the coincidence limit for the arguments, from the corresponding formula for the 
Wightman function we obtain the boundary-induced part in the VEV of the field square: 



<P0 ^ Jm K qn {za) 
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0.5 1 1.5 2 2.5 0.5 1 1.5 2 



Figure 3: Boundary-induced parts of the components of the energy-momentum tensor multiplied 
by a D+1 , a D+1 (T-) a , as functions of r/a in both interior and exterior regions of the cylindrical 
surface. The graphs are plotted for minimally (left panel) and conformally (right panel) coupled 
D = 3 massless scalar fields with Dirichlet boundary condition on the cylindrical surface and 
for the cosmic string bulk with q = 4. The numbers near the curves correspond to the value of 
the index i. 



with Ad defined by Eq. (|4.4|) . As for the interior region, the expression on the right diverges 
on the cylindrical surface. The leading term in the corresponding asymptotic expansion near 
this surface is obtained from that for the interior region, formula (|4.6|) . replacing (a — r) by 
(r — a). For large distances from the cylindrical surface, r S> a, and for a massless scalar field, 
we introduce in 1)5.1 J) a new integration variable y = zr/a and expand the integrand over a/r. 
The main contribution comes from the n = term. By taking into account the value for the 
standard integral involving the square of the MacDonald function 02j, to the leading order for 
A ^ 0, one finds 

n l-D/2 p2 (D^l) 

^ )a ~- 2 D^r D -Hn{r/a) r(D/2) ' (5-2) 

In the case of Neumann boundary condition {A = 0) and q > 1, the leading contribution again 
comes from the n = term: 

IA . ^!3i (-r. 



ja^- 1 r(D/2) V2r 

For Neumann boundary condition and q < 1, at large distances the main contribution comes 
from the n = 1 term and the VEV of the field square behaves as l/r^ 2 * 3-1 . As we see, the 
boundary-free part in the VEV dominates at large distances from the boundary. For a massive 
field, under the condition mr 3> 1 the main contribution into the integral in 1)5.1)1 comes from 
the lower limit and one obtains 



oo 



(4vr) — o r^- ^ K qn (ma) 

with the exponential suppression of the boundary induced VEV. Note that here the suppression 
is stronger compared with the boundary- free part. For large values 5>1, the contribution of the 
terms with n ^ to the VEV 1)5.1)1 can be estimated by using the uniform asymptotic expansions 
for the modified Bessel functions when the order is large. This contribution is suppressed by 
the factor q( D ~ 1 ^ 2 (a/r) 2q and the main contribution to the VEV of the field square comes from 
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the n = term which is a linear function on q. In the opposite limit, q <C 1, analogously to 
the procedure for the interior region, it can be seen that {(p 2 ) a tends to the finite value which is 
obtained from (|5.1jl replacing the summation over n by the integration. 

For the part in the vacuum energy-momentum tensor induced by the cylindrical surface in 
the region r > a, from (|2.22[) . ()3.14|) . (|5.1[) one has the following formula 

^ = ¥ E'T dzZ " & - ™ 2 ) ^ TFT^S i K ^M] , (5.5) 
00 ^ i m K qn (za) 

with the functions F$ [f(y)\ defined by formulae (jUnj-glSI). The VEVs given by formula 
(|5.5jl diverge at the points on the bounding surface with the leading divergence obtained from 
analog formula in the interior region, Eq. (|4.21j) . by the replacement a — r — > r — a. By the way 
similar to that used above for the vacuum expectation value of the field square, it can be seen 
that for a massless scalar field with A ^ and at large distances from the cylindrical surface, 
r 3> a, the components of the vacuum energy-momentum tensor behave as 

(T °) a « -(D - l)(Tl) a « ) a ~ (47r) D/2-i 2 r( ^ /2) 0or i+i ^ r/a) ' ( 5 - 6 ) 

For a conformally coupled scalar field this leading terms vanish and it is necessary to take 
into account the next terms in the asymptotic expansion. In the case of Neumann boundary 
condition the components of the energy- momentum tensor behave as l/r D+3 when q > 1 and as 
\/r D+2( i +1 when q < 1. For a massive scalar field the main contribution comes from the lower 
limit of the integral in Eq. (|5.5[) and we find 



2 i?-l 7r (l?-3)/2 V r / ^ K 9 „(ma 

As we see in this case the radial stress is suppressed compared with the other components by an 
additional factor mr. As in the case of the field square, for large values of the parameter q S> 1 
the main contribution to the VEV of the energy-momentum tensor comes from the n = term 
which is a linear function on q. For small values q <C 1, in the leading order the summation 
over n can be replaced by the integration and the VEV tends to the finite limiting value. In 
figure we have plotted the dependence of the boundary-induced parts in the components of 
the vacuum energy- momentum tensor for the exterior region as functions on r/a in the case of 
D = 3 Dirichlet massless scalar with minimal and conformal couplings. 

In the discussion above we have considered the idealized geometry of a cosmic string with 
zero thickness. A realistic cosmic string has a structure on a length scale defined by the phase 
transition at which it is formed. As it has been shown in Refs. |141 123j . for a non-minimally 
coupled scalar field the internal structure of the string has non-negligible effects even at large 
distances. Note that when the cylindrical boundary is present with the boundary condition 
(j2.3|) . the VEVs of the physical quantities in the exterior region are uniquely defined by the 
boundary conditions and the bulk geometry. This means that if we consider a non-trivial core 
model with finite thickness b < a and with the line element (|2.1[) in the region r > b, the results 
in the region outside the cylindrical shell will not be changed. As regards to the interior region, 
the formulae given above are the first stage of the evaluation of the VEVs and other effects could 
be present in a realistic cosmic string. 

6 VEVs for a cosmic string with finite thickness 

From the point of view of the physics in the exterior region the cylindrical surface with boundary 
condition (|2.3j) can be considered as a simple model of cosmic string core. In general, the string 
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core is modelled by a cylindrically symmetric potential whose support lies in r ^ a. In this 
section we generalize the results in the exterior region for the model of core described by the 
line-element [121123] 



A? 



ds 2 



dt 2 



P 2 (r/a)dr 2 - r 2 



(6.1) 



i=l 



where P(x) is a smooth monotonic function satisfying the conditions 

lim P(x) = 1/q, P{x) = 1 for x > 1. 
x— »0 



3.2) 



The eigenfunctions in the region r < a have the structure (|2,5|) with the radial function 
f n (r/a,^/a) instead of P a Jq\n\(l r )- The equation for the radial function is obtained from field 
equation ()2.2[) with the metric given by ()6.1|) : 



d x d 22 q 2 n 2 
+ 7 a 



xP{x) dx P{x) dx 



2£ P'(x) 
x P 3 (x) 



/ n (x,7o) = 0. 



(6.3) 



We will denote by R n (x, 7a) the solution of this equation regular at r = 0. Now the radial part 
of the eigenfunctions is written in the form 



R n (r/a,~/a) for r < a 

A n J q \n\ il r ) + B nY q \ n \ (jr) for r> a. 



(6.4) 



The coefficients A n and B n are determined from the conditions of the continuity of the radial 
functions and their derivatives at r = a: 



Ar 



TT 

2 



TV 

2 



7«^'|„|(7«)^n(l,7«) - Y q \n\(ja)R' n {l,>ya) 
7aJg| n |(7«)^n(l,7«) - J g \n\(ia)R' n (l,-fa) 



(6.5) 
(6.6) 



where i?^(l,7a) = (d/dx)R n (x, ^a)\ x= \. Substituting these expressions into the formula for the 
radial eigenfunctions in the region r > a, we see that these eigenfunctions have the form (|2.17|) 
where the barred notations are obtained from from the expressions given by Eq. (|2.8f) by the 
replacement A/B — > —R' n (l,^a)/R n (l,ja). The further evaluation of the Wightman function 
in the region r > a is similar to that described in subsection 12.31 Therefore, the part in the 
Wightman function induced by the non-trivial structure of the string core is given by formula 
()2.22() . where in the expressions for the definition of the barred modified Bessel functions [see 
Eq. (|2.<S|) ] we should substitute 

A_ R' n (l,zae^ 2 ) 

B R n (l,zae™/ 2 )' 1 ' ' 

The formulae for the VEVs of the field square and the energy-momentum tensor are obtained 
from formulae (|5.1|) and (|5.5|) by the same substitution. The corresponding results for the 
Euclidean Green function and the VEV of the field square in a special case of D = 3 massless 
scalar field are given in Ref. |23J. Two specific models of the string core have been considered in 
literature. In the "ballpoint pen" model the region r < a has a constant curvature, whereas 
in the "flower pot" model |14j the curvature of spacetime is concentrated on a ring of radius 
a. In these models the Euclidean Green function and the VEV of the field square for D = 3 
massless scalar field are investigated in Ref. |14| . For the first model one has 



1 



Pl n \l/q)' 



-2£ 



z 2 a 2 



(6. 
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where P„ (x) is the associated Legendre function. For the second model the curvature is a 
delta function concentrated on a ring of radius a and the radial parts of the eigenfunctions 
have a discontinuity in their slope at r = a. The corresponding jump condition is obtained by 
integrating the radial part of the field equation through the point r = a. This procedure leads 
to the ratio of the coefficients given by the formula 

A I[ \(za/q) 

^ = -za-^-^-2t( q -l). (6.9) 
B I\ n \{za/q) 

For the first model in (2+l)-dimensions the ground state energy of a massive scalar field is 
investigated in Ref. 



7 Conclusion 

We have investigated the local one-loop quantum effects for a massive scalar field induced by a 
cylindrical boundary in the spacetime of a cosmic string. We have assumed that on the bounding 
surface the field obeys Robin boundary condition. The latter is a generalization of Dirichlet and 
Neumann boundary conditions and arises in a variety of physical situations. As a first step in 
the evaluation of the renormalized VEVs of the field square and the energy-momentum tensor, 
in section [21 we have considered the Wightman function in both interior and exterior regions. 
The corresponding mode-sum in the interior region contains the summation over the zeros of a 
combination of the Bessel function and its derivative. For the summation of the corresponding 
series we have used the generalized Abel-Plana formula which allows us to extract from the 
mode sum the Wightman function for the cosmic string background without the cylindrical 
shell and to present the boundary-induced part in terms of exponentially convergent integrals 
in the coincidence limit for the points away from the boundary. The representation of the 
Wightman function where the boundary-free part is explicitly extracted is given also for the 
exterior region. The boundary-induced parts in the interior and exterior Wightman functions 
are related by the replacements I qn T± K qn . 

In section 01 we have considered the VEVs induced by the cosmic string geometry without 
boundaries. Though this geometry is well-investigated in literature, to our knowledge, no closed 
formulae were given for the VEVs of the field square and the energy-momentum tensor for 
a massive field with general curvature coupling in an arbitrary number of dimensions. We 
show that such formulae can be derived when the parameter q is an integer. In this case the 
corresponding Wightman function is the image sum of the Minkowskian Wightman functions. 
Renormalized VEVs of the field square and the energy-momentum tensor are determined by 
formulae (|3.7|) . (|3.15f) - (j3. 17f) for a massive field and by formulae (|3.1(Jj) . 1)3.22(1 for a massless 
one. In the latter case and for odd values of the spatial dimension the summation over I can be 
done by using the recurrent formula. In this case the VEVs are polynomial functions of q and 
by analytic continuation the corresponding formulae are valid for all values of this parameter. 
By using the formula for the interior Wightman function, in section 0] we have investigated the 
VEVs of the field square and the energy-momentum tensor in this region. The corresponding 
boundary-induced parts are given by formulae Q4.3JI and l|4.12|) . For the points on the string 
these parts are finite when q ^ 1. For q < 1 the boundary- induced part in the VEV of the 
field square remains finite on the string, but the corresponding part in the energy-momentum 
tensor has integrable divergences. Near the string the boundary-free parts behave as \jr D ~ x 
for the field square and as l/r D+1 for the energy- momentum tensor and these parts dominate. 
For the points near the boundary the situation is opposite and the boundary-induced parts are 
dominant. For large values of the parameter q the boundary-induced VEVs for both field square 
and the energy-momentum tensor are linear functions of this parameter, whereas for small values 
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of q they tend to finite limiting value. We have the similar behavior for the VEVs in the region 
outside the cylindrical shell. These VEVs are investigated in section |S] and are given by formulae 
(|5.1jl and (|5.5|) . In the case of a massless field with non-Neumann boundary condition, at large 
distances from the cylindrical surface the boundary-induced parts behave as 1/ [r D_1 ln(r/a)] 
for the field square and as 1/ [r B+1 ln(r/a)] for the energy-momentum tensor. For a massive 
field under the assumption mr ^> 1, the boundary induced VEVs are exponentially suppressed. 
This suppression is stronger than that for the boundary-free parts. 

The cylindrical surface with boundary condition (|2,3j) can be considered as a simple model of 
the cosmic string core. In section|H]we give the generalization of the corresponding results in the 
exterior region for a general cylindrically symmetric static model of the string core with finite 
support. We have shown that the corresponding formulae are obtained from the formulae for 
the cylindrical surface with the substitution given by Eq. I|6.7|) . where R n {r /0,7a) is a regular 
solution of the radial equation for the eigenfunctions. For two special models of the string core, 
namely, the "ballpoint pen" and "flower pot" models, the ratio of the coefficients is given by 
formulae (|6.8|) and (|6.9|) . respectively. 
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